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and classified. The Lagrange multipliers are derived. It is shown that a 
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in general. Both modes exhibit a simple type of "constraint bifurcation" for 
certain special field/parameter values. 
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1 Introduction 

The Poincare gauge theory of gravity (PGT)lila, based on a Riemann-Cartan ge- 
ometry, allows for dynamic torsion in addition to curvature. Because of its gauge 
structure and geometric properties it was regarded as an attractive alternative to 
general relativity (GR). It was quickly realized that the theory had physical difficul- 
ties with generic values for its ten coupling parameters (see, e.g., referencesoU'ofi'Q) . 
Consequently investigators have looked for restrictions giving viable sets of PGT pa- 
rameters. By utilizing certain theoretical tests (e.g., "no-ghosts" or "no-tachyons"), 
sets of constraints on the parameters and possible viable PGT modes were obtained 
(see, e.g., referencesBH'Qla'Ll) . These investigations naturally used the weak- field ap- 
proximation and linearization of the theory, avoiding the inherently highly nonlinear 
complications in the PGT. 

Subsequently, investigations^ indicated that certain degeneracies of the Hessian 
matrix, thought to be necessary for a viable theory, appeared to lead to difficulties 
with the initial value problem. Then a shock-wave analysis^ty concluded that these 
same degeneracies allowed tachyonic propagating modes. These difficulties were 
more recently reconsideredtll. It was found that, when all the constraints are taken 
into account, there are no problems — at the linear order. 

However, the PGT Hamiltonian analysis tells a more complicated story. The 
virtue of the Hamiltonian analysis is that it provides a clear vision of the possible 
degeneracies, corresponding constraints, and true degrees of freedom of a theory. 
Its application to the nonlinear PGT is revealing. It shows that the nonlinear 
behavior of the PGT can be — and, through a phenomenon referred to as "constraint 
bifurcation" t3, very likely will be — qualitatively different from the linearized one in 
the number and type of constraints, so the linearized "good modes" may very well 
not be viable in the full nonlinear theoryll3. Hence, in order to understand the subtle 
behavior of the PGT and search for modes which truly have good propagation, it is 
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important that the analysis consider the full nonlinear scope of the theory. 

It has long been known that non-linear couplings involving higher spins are 
problematicalO'Ej. However, spin-zero modes were expected to be problem free. For 
the PGT, it has subsequently been verified that two four-parameter subclasses with 
only spin-zero propagating modes really do have a well posed initial value problem 
without any tachyonic propagationll3. 

That leads us to examine the Hamiltonian formalism for some representative 
dynamic spin- zero PGT modes in order to see how (and indeed whether) the modes 
really manage to avoid the nearly ubiquitous (and almost certainly fatal) nonlinear 
"constraint bifurcation" . 

Such an investigation seems to be an essential prelude to a search for any higher 
spin nonlinear-problem-free modes in the PGT. Our results presented here gives, to 
our knowledge, the only complete Hamiltonian analysis of PGT modes believed to 
be (linearly and nonlinearly) viable. 

Also, it may be worthwhile to mention, certain interesting results which were 
found in studies of higher-derivative gravity (see, e.g., referencescHMEJlHy) . These 
theories contain a non-ghost scalar field, in addition to the usual graviton of GR. 
Some pure gravity inflationary models for the Universe were proposed based on such 
theories, although the issue still remains controversialEil'E3Ho. This has provided 
additional motivation for looking for similar situations in the PGT, but it must be 
kept in mind that the principle and structure between the PGT and higher-derivative 
gravity are quite different. 

The paper is organized as follows. In section 2 we review the basic elements of 
the PGT and introduce its Lagrangian and field equations. In section 3 we use the 
Dirac theory for constrained Hamiltonian systems in the excellent "if" constraint 
formulation developed by Blagojevic and NikolicB§B3. The primary constraints, in- 
cluding ten "sure" primary constraints and thirty so-called "primary if-constraints" , 
are found. The total Hamiltonian density, including the canonical Hamiltonian den- 
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sity and all possible primary constraints, is derived. In section 4 we consider two 
very degenerate spin-zero modes, each having only one of the six parameters of the 
quadratic curvature parts being non-zero. In particular, we consider a spin-O"*", a 
massive spin-0~ and the associated massless spin-0~ modes. The Lagrange multi- 
pliers are derived. It is shown that the massless spin-0~ mode is unphysical, being 
pure gauge. In section 5 the degrees of freedom are counted. The multipliers are 
shown to be exactly the missing 'velocities'. Their effects in the spin- zero modes 
are discussed. In both scalar modes we find a simple type of "constraint bifurca- 
tion" phenomenon, which is connected with changes in the nature of the constraint 
reduced from the Lorentz rotation parts of the canonical Hamiltonian density. The 
similarity between the spin-0"'" case and higher-derivative gravity is noted. In the 
final section we present our conclusions. 

Throughout the paper our PGT conventions are basically the same as Hehl'sll. 
We have made a few adjustments to accommodate the translation of the Hamiltonian 
"if" constraint formalism to these conventions. The latin indices are coordinate 
(holonomic) indices, whereas the greek indices are orthonormal frame (anholonomic) 
indices. The first letters of both alphabets (a, b,c, . . .; a, /3, 7, . . .) run over 1, 2, 3, 
whereas the later ones run over 0, 1, 2, 3. Furthermore, r]^^ =diag(— , -|-, +, +); e'^^'^^ 
is the completely antisymmetric tensor with e°^^^ = — 1. The meaning of a bar over 
a greek index is adopted from Blagojecvic and Nickolic c3. 

2 Poincare gauge theory of gravitation 

In the PGT there are two sets of gauge potentials, the orthonormal frame field 
(tetrads) et'^ and the metric-compatible connection Fi^'^, which are associated with 
the translation and the Lorentz subgroups of the Poincare gauge group, respectively. 
The associated field strengths are the torsion 

T,/ = 2(9[,e,/ + r[,|/e|,/), (2.1) 



and the curvature 

R^j^J.'' = 2(<9[jrj]^'' + r[i\^T\j]^") , (2.2) 

which satisfy the Bianchi identities 

"^[iTjkf = R[ijkf, (2.3) 

V[^R,k^ = 0- (2.4) 

The conventional form of the action, which is invariant under the Poincare gauge 
group, has the form 

A = Jd^xeiLM + Lc), (2.5) 

where Lm stands for the matter Lagrangian density (which determines the energy- 
momentum and spin source currents), Lq denotes the gravitational Lagrangian den- 
sity, and e = det{ei^). In this paper we are concerned with the gravitational propa- 
gating modes, hence we omit the matter Lagrangian density, so Lq is considered as 
the source-free total Lagrangian. Varying with respect to the potentials then gives 
the (vacuum) field equations, 

V,V;' - e; = 0, (2.6) 

V,P^.^^' - e^,' = 0, (2.7) 



with the field momenta 



and 



V^^ ■= M^ = 2^^ (2 8) 

,, ._ deLc _ deLc 



e; := e\eLG-T^,''VJ'-R^r'P.J\ (2.10) 

e^J ■■= V[,^]\ (2.11) 
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The Lagrangian is chosen (as usual) to be at most of quadratic order in the field 

strengths, then the field momenta are linear in the field strengths: 

e ^ C") 
V = 7iE«^^^. (2-12) 

V,.'^ = -^e\,e\^ + -J2b,W\,, (2.13) 

' '* fc=i 

W.. (fc).. 

the three T-'*^ and the six W^ ^^ are the algebraically irreducible parts of the torsion 

and the curvature, respectively. The reciprocal frames e*^ and Cj^ satisfy e^ci" = 5^^ 

and e\ej^ = 6/; the coordinate metric is defined by gij = ei^Cj^rj^y. The a^ and 

hk are free coupling parameters. Due to the Bach-Lanczos identity only five of 

the six 6fc's are independent, a^ is the coupling parameter of the scalar curvature 

R := R^u^^. For the Hamiltonian formulation we associate the canonical momenta 

with certain components of the covariant field momenta: 

vrV = V;\ (2.14) 

ttV. = ^M-*. (2.15) 

3 Primary constraints and total Hamiltonian 

In this section we present the primary constraints and the total Hamiltonian density 
of the PGT in terms of the decomposition of the canonical variables and fields. First 
of all, one obtains the "sure" primary constraints 

A ~ 0, (3.1) 

T^\. ~ 0. (3.2) 

These constraints just refiect the fact that the torsion and the curvature are defined 
as the antisymmetric derivatives of Cj^ and Fi^^; they do not involve the "velocities" 
Cq^ and Fq^^ . One will obtain further primary constraints if the Lagrangian density is 
singular with respect to the remaining "velocities" , Ca^, and Fa^^ . (Such constraints. 
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so-called "primary if-constraints" , result from certain vanishing coupling parameter 
combinations.) The total Hamiltonian density is of the form 

^ot = ^can + Uo^TT^ + \^urT^\u + u^<\>A, (3.3) 

where the 0^ are the primary "if" constraints and the m's denote the associated 
Lagrange multipliers; Tican stands for the canonical Hamiltonian density which will 
be specified below. 

Before we proceed to obtain the explicit form of the canonical Hamiltonian den- 
sity and (pA, it is convenient to define the decomposition of related variables and func- 
tions. We essentially follow the techniques developed by Blagojevic and NikolicE3. 
Let us note that the components of the unit normal n to the x^ =constant hyper- 
surface, with respect to the orthonormal frame, are given by 

'V - ^. (3.4) 

A vector, e.g., V^, can be decomposed naturally into the orthogonal and parallel 
components with respect to the orthonormal frame indices: 

V^ = -V^n^ + Vp, (3.5) 

Vi. = VX, (3.6) 

V-j, = K(V + ^X)- (3-7) 

One can easily extend the decomposition to any tensors with orthonormal frame 
indices. The lapse and shift functions can be written as 



1 



N = ^^ = -n^eo^ (3.8) 



Oa 

N^ ^ -^ = eoV, (3.9) 

and e = NJ, where J is the determinant of the 3-metric. Defining the convenient 
"parallel" canonical momenta 

Tt'^^ = Ca'^TT^fj., n^'fj.u = ea'^Tr^fj.u, (3.10) 



which satisfy tt'^ ^Ua- = 0, tx'^ ^uf^a = 0, the canonical Hamiltonian density, 

Hcan = TT^e/ + ]^'K\,rr - eL, (3.11) 

can be rewritten in the so-called Dirac-ADM 

Hcan = Nn± + N'^Tia + \rrn^u + daV\ (3.12) 

which is linear in A^ and N"-. The other quantities are given by 

H± = 7r%Tx^'^ + ^TT Vi^x^'^'^ -JL- n^Va7T\, (3. 13) 

Ha = TT^Tab^ + ^TT^.i^ab'^" ' ea^V^TT^, (3.14) 

Only the super-Hamiltonian 7i± is involved in dynamical evolution, the super- 
momenta Tia and Lorentz rotation parts 7^^,^ are kinematical generators, conse- 
quently we concentrate on 7i^ when consistency conditions are calculated. By uti- 
lizing the forms of the total and canonical Hamiltonian density ( |3.3|) and (|3.12|) in 
the consistency conditions for the primary constraints ( p.lj ) and ( pj.2| ) we obtain the 
secondary constraints (SC) 

^±-0, Ha^O, n^u^O. (3.17) 

Since it is easy to check that the primary constraints ( p.l| ) and ( |3.2| ) are first-class, 
i.e., Tc^fj, and vr^^i, are unphysical variables, by using the Hamilton equation of motion 
we can infer that the multipliers Uq'^ and Uq'^^ are indeed equal to Cq^ and io^'^. They 
are dynamically undetermined pure gauge multipliers. 

It is necessary to understand the relation between the canonical momenta and 
the velocities well before one can make 7i± more apparent. Let us consider the 
torsion momenta ir-pi, at first, n-pu can be decomposed into four irreducible parts as 



follows: 

= -n^'n^-pj_ + 7T-pu + 7T-pjj+ -rjjjjrTl, (3.18) 

A '■^ 

where vr^, it and tt^ are the antisymmetric part, trace part and symmetric-traceless 



part of TTjjij, respectively. Manipulating the definition of the torsion momenta ( p.l4|) 



one finds the following relations between the different parts of the canonical momenta 
and the corresponding parts of the velocities T±-p^: 






V± = —r- ^7?(ai - (^2)TjT = -^7^(2ai + a2)T^jz±, (3.19a) 



A n- 1 



'W- —r~ ^("1 ~ a3)Tjiu± = -^(«i + 2a3)r^[7ij7], (3.19b) 



J 3/2' ^ '^ 3/ 



Ml' T /2 



J /2 

TT _ 02 



T±{jn7), (3.19c) 



TfT^a", (3.19d) 



where T^ = Tj^j^, and a tensor with two indices contained in the bracket ( ) de- 
notes that the tensor is symmetric-traceless with respect to the two indices. If the 
parameters take on the critical values: 2ai + 02 = 0, ai -|- 203 = 0, ai = 0, and/or 

A 

02 = 0, one obtains the following ''primary if- constraints" (PIC): 0^_l ^ 0, 0^ ~ 0, 
^jiTT ~ 0) and/or 0^0, respectively. 

The curvature momenta TTa^u can be decomposed into six irreducible parts: 

^!'a^lu = T'wjTiJ + 27r^±i-pnu], (3.20) 

and 

A '~^ -L / \ 

T^-pUL = T^JTu^ + ^A«7± + -^^jnjT^ ±-, (3.21) 

1 __. 4^ 

T^-dw = --e<T^!.±% + T^ijirjup + -''tts^Imf], (3.22) 

where the notation is as follows: For a spatial tensor Xg^ = Xg:[^], ^X = e'^'^'^^X^^, 
Xj[ = Xtj/, and X^^ = X(;^\-p\j,) - \X-pijj:^ + ^rfj^pX^) are the pseudoscalar part, 
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vector part, and traceless tensor part of Xg^, respectively. Identifying the irre- 



ducible parts of the curvature momenta (|2.15D, one finds 



a^iv 



T + -(^2 - hYRo± = --(&2 + hYR±o, (3.23a) 

J K K 

^ - -(&4 - h)Rj,^ = -(64 + h)R^j,, (3.23b) 

^ - -(&1 - h^YRjj^^ = --(61 + h^YRi.^, (3.23c) 

J K K 

Cj)^= ^ + ^ + T^(&4 - h)R = -\W + &6)i?±±, (3.23d) 

(Pw^ - ~^ + ~(^2 - h)R[jip] = --(^2 + b5)R±[jiu]±, (3.23e) 

+ -(&i - h)Rg^) = -i(6i + 64)i?±(7iF>±, (3.23f) 



W = ^p ^ K*-^^ ~ bi)R{iiu) ^ 



where Pi?o± := e'^'^'^Rj^^, m^o := e'^""^^?^^^, E^zf := ^^", and E := E/- By a 
similar argument as used above, for various degenerate parameter combinations one 
can obtain any of the six expressions of (|3.23a| -f) as PIC's. The relations between 



the critical parameter combinations and the constraints are summarized in table 1. 
In order to treat all such possibilities in a concise way, the singular function 

\(x) I l/x, X ^ 0, 

^^^ = <^ ^ ^ (3.24) 

^ [0, X = 0, 

was introduced. The PIC's in the total Hamiltonian density ( |3.3| ) can be given in 
the form 

u^^A = {u ■ 4>f + {u ■ 0)^, (3.25) 

where 

{u ■ (t)f = [1 - A(2ai + a2)]M^^07i± + [1 - \{a^)]u^^^ 



A— ^ 1, 

'""6—4- 

3' 



+ [1 - A(ai + 2a3)]uf"'(t)j^ + ^[^ " H(^2)W (3-25a) 



and 



{U ■ 0)^ = ^[1 - \{b2 + h)ruP<P + ^[1 - Kh + 62)]^«'^^^^0^ 
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+[1 - X{h + h)]u^<Pj, + 2[1 - X{b2 + h)]u^''^<jyj^^ 
+2[1 - A(6i + b,)]u'^^4>w± + ^[1 - ^(^4 + &6)]w^0±. (3.25b) 
The super-Hamiltonian 7i± in Tican ( P-12 ) then turns out to be of the form 



7Y_L = 7i\ + 7i 



±) 



(3.26) 



with 



where 



n] 



n 



R 



-Jl' 



3A(2ai + 02) 
2ai + 02 



>fi±^ 



3A(ai + 203) A A, 



ai + 2a3 



'TZ^'^^ 



4 

4 



'^(«i);_:77F^ A(a2)^2 



|1U^ 



tti '" 3a2 

A(&2 + h) „^2 , ^(64 + h 



Jl7 - n^Va7T\, 



-PA,- + 



24(62 + &3)" ' 4(64 + 65)"'^ 



3(61 + 62) ""'^^ 

A(6i + M7 ^ 



+ 



A(62 + &5) ;;)_ ^ 

2(62 + ^5)'^^'^^ 



2(61 + 64) 



^77i7±S 



i/xl/± 



+ 



A(&4 + K) 
6(64 + &6) 



^±<( 



(3.26a) 



- JL!", (3.26b) 



ij 



12/2 



(2ai + a3)Tj;^^T'^'^'^ 
+2(ai - a3)T^r^ - 2(ai - 02)?^?^ 






^/lU- 



(3.26c) 



A is the cosmological constant. The relations between the constants Cj's 
6j's are given by 

1 



Co 

Cl 
C2 



-^(261 + 362 + 63), 
-7^(&i-&3), 



6k 
1 



(261-362 + 63), 



(3.26d) 
and the 

(3.27a) 
(3.27b) 
(3.27c) 
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C3= ^ibi + b2-h-b,), (3.27d) 

C4= ^(6i-&2-&4 + 65), (3.27e) 

C5= -77^(261- 364 + &6). (3.27f) 

We will now apply this wonderful general "if-constraint" Hamiltonian formulation 
to certain spin-zero modes of PGT with specific parameter combinations. 



4 Spin-zero modes 

The PGT propagating modes with a single spin-zero propagating mode are sup- 
posed to have ghost-free and tachyon-free charactergHlfl'Qla'Ea. There are two spin-zero 
modes, i.e., 0"*" and 0^ according to table 1. These spin- zero modes were analyzed 
by a covariant Lagrangian technique and were shown to have a well posed initial 
value problem with no tachyonic propagation characteristicsts. In order to better 
understand them, and the non-linear PGT problems, it is worthwhile to examine 
their full nonlinear behavior under the Hamiltonian analysis and how (and indeed 
whether) they really avoid the "constraint bifurcation" problems (which could be 
expected to arise from the many nonlinear constraints). 

The restrictions ai + 2a3 = and 2ai + 02 = have been regarded as "viable" 
conditions for the PGT theories, we will likewise (and in the interests of simplicity) 
also assume these restrictions in both of our propagating spin-zero parameter choices. 

4.1 Spin-0^ mode 

According to Table 1. 7t± corresponds to spin-0^. We make the specific parameter 
choices: 

2ai + 02 = 0, ai + 2a3 = 0, 

hi = h2 = h = h = h^ = 0, (4.1) 
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ao ^0, foe ^ 0. 



In fact it is not necessary to take 6i, . . ., 65 to vanish. The choice can be relaxed 
by requiring only 61 = —62 = &3 = — &4 = &5 7^ (refer to ( |3.23a[ -f)). However, we 
believe that it's more critical to understand the character of the pure spin-0'*' mode. 
This simple choice will greatly simplify calculations, while more general choices are 
expected to have the same qualitative behavior. After all the 4-covariant analysisllj, 
has argued that there is a 4-parameter class of Lagrangians with a dynamic spin-0"'", 
but each has qualitatively the same dynamic behavior. Hence, we expect that the 
Hamiltonian analysis of these more general choices, aside from being calculationally 
much more complicated and less transparent, would actually show no new interesting 
dynamic features, so we leave them for future work. 
The corresponding super-Hamiltonian is: 



H1 



T^T + ^f +, 



(4.2) 



and 



n' 



H 



R-\ 



l" 



2Jai 



vr^vr^'^ + 



l2Jai 



-n^ - n^VcT^\ - JU 



^.tLl 3ao _ _^ox2 
666 J l^ 2k-^ 

1 rp rpT7al_i I ^1 r p rpfliju 



||tf.|^.«-.A, 



'^1 r T^ 
2P '^ ■ 



(4.2a) 
(4.2b) 
(4.2c) 



Due to the parameter choice ( ^.1| ), the PIC's are as follows: 



A 



A 

J 2P ''"^ ' 

% 
-«0, 



^ ~0, 



J 



(4.3) 

(4.4) 
(4.5) 

(4.6) 



'(t>WjM7 = 


J 


^0, 


A 


A 
J 


^0, 


(t)-rTT7\ = 


'^]I17± 


^0. 
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(4.7) 

(4.8) 

, ^. (4.9) 

According to Dirac-Bergmann algorithm it's necessary to identify the class of these 
constraints. The non-vanishing Poisson brackets (PB) for the PlC's are the follow- 
ing: 

A A (5 , 

{</'77F±,0W} ~ -^^V^[-pVu]r, (4.10) 

{<PjT,<Pv±}- — ^^rfj^- (4.11) 

A 



where w := —^ -\ — --. Thus (f)ji±, (f)-pjj, 0^ and 0jjj7_|_ are second-class — as long 



7r_|_ 3ai 
J P 
as ro 7^ 0. The constraints ^0, '^(fy^jju and (f)jxu± commute with the other primary 

constraints; according to the general PGT "if" constraint analysisE3'E3, the associated 

SC constraints can be derived from them by calculating their time derivatives. They 



are 



P0 ^ --Nw^T ^ 0, (4.12) 

^\nw^^O, (4.13) 



^Ja7± 



s. 



3ai J 

'Wjn7^ —Nw^j^^O. (4.14) 

For the time being we put aside the highly degenerate case tu ^ 0. Then, assuming 
that m does not vanish (except perhaps on a set of measure zero), these SC's can 
be simplified to 

Px= T ^ 0, (4.15) 

Xjiu±^ ^^0, (4.16) 

'Xajn7^ T^^O. (4.17) 
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By identifying and eliminating all primary and secondary constraints inside, the 
supermomenta and Lorentz rotation parts lead to the distinct constraints: 

K^Q => -yi?a±-5„^ + ^TSV6e/^0, (4.18) 

n^^^O => wTj^^^O, (4.19) 

ru % +e%dazu ^ 0. (4.20) 

All of the "if" constraints (including the PIC's and the SC's) are eventually 
second-class by calculating their PB's with the others. The details are shown in the 
appendix. The consistency conditions have to be obeyed, i.e., we must produce the 
time derivatives of second-class constraints and force the results to vanish weakly. 
We then learn the Lagrange multipliers. The procedure can be understood from the 
following expression: 



'^B= J {(pB, Hl^^} 

= /[{0B,Kan} + ^'^{0i?,</>A}] ^0. (4.21) 



If (f)B is second-class, then it is possible to solve for some u . In Titot (|3.3| ) only 
the super-Hamiltonian 7i± (namely 7i|[| here) involving time evolution is important. 
By proceeding with the step described by ( [4.21| ), it is straightforward to get the 
Lagrange multipliers: 



u>' 



u^"^ = --NJT^", u^"" = 0, (4.22a) 

-NJR^^, u'^^ = --NJBp^\ (4.22b) 

% = ^Njm^u (4.22c) 

%^^ = -NrW^^, (4.22d) 



u^""^ = -NJ 






(4.22e) 
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Now the canonical Hamilton equations of motion can be derived directly from 
the completed total Hamiltonian density, 



qA= JiqAMot}, (4.23) 

^A= J{7rA,n'tot}, (4.24) 

where q^ represents the collection of canonical variable Cj^ and Fi^'^, tia represents 
the collection of the conjugate momenta vr*^ and Tr*^,^. 

4.2 Massive spin-0~ mode 

From table 1, % corresponds to the spin-0^ dynamic mode. We consider the simple 
parameter choice: 

2ai + 02 = 0, Oi + 2a3 = 0, 

b,=b2 = h = h = be = 0, (4.25) 

ao ^0, 63 7^ 0. 

The corresponding super-Hamiltonian is: 

ni= nl + n^-, (4.26) 



and 



24b':!. J K 24k 2r 



nr = -irri^ - -"Roi.? - ^m^^m^ + J,jr - ja. (4.26a) 



In this case the PIC's are: 



J P 



<h^^ ^-^Tj,^ 0, (4.27) 



A 
TTtju Ol 



W= — - ^^A-± ^ 0' (4-28) 
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^ ^ 0, (4.29) 



J 



' a^u 



A 



J 

7r_L 3ao 
J / 



0, (4.30) 



^ + ^ - 0, (4.31) 



A 

^ ^ 0, (4.32) 

^^^ ^ 0. (4.33) 



J 
The non-zero PB's for the PIC's are the following: 

A 5^^, % 

{(t>JIu±^(t>'ai}- -^-j^^va±, (4.34) 

A A ^ , 

{'/'tjj^X' 0W} ~ — ^"^^[ji^]-, (4.35) 

{0M,0^±}^ 2^mr^, (4.36) 

where m = — -5 — . We first consider the generic "massive" case: m 7^ 0. Since the 
determinant of the PB's matrix is J^"'^^[(%/6J)^ + m^]^ > 0, the constraints 0jj_l, 

A ^ A 

0^, 0^ and (f)ju;^ are second-class. 



Ju7± aLL>u (pcrfiu^ 



Similarly we have the time derivatives of 0_l, (p-pjj^ and 

: 1 _% /^ TT 



-- l2^7+'"a;7'=° <"«' 



/2 % 



^Mi7± - Q T^C^r'^^^'^l^) + ^-^ ^ 0' (4-39) 

y J cLi J 



"(p-^jiu ~ — ^e^M(^|_L7r|5.)^ + mTjjj^ ^ 0. (4.40) 

Generically these can be replaced by the three simpler SC constraints: 

1 % f -K 

X±= -r7;'T- + m—-^0, (4.41) 

iz ,j di J 

Xj^i.^ ^ ^ 0, (4.42) 

\wfn7 = T^ ^ 0. (4.43) 
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The supermomenta and Lorentz rotation parts are: 



'^1 rr!vm__ _|_ '^1 pji rpvW 



"-L rjiyrj i 



1 % 
127 



^tio ,n -'^/y. 



1-/10 



/2 



l-JI± 



/2 



T^'Van^ 



TT 



^M^^O 



1 % — ^ 

1 _ % %^- 



(4.44) 
(4.45) 
(4.46) 



where ^Rjio := e'^'^'^^RjuPfw- Utihzing again the standard Dirac-Bergmann algorithm 
the corresponding Lagrange muhiphers are found to be 



u 



-pL 



W 



u 



s^cyi^v 



U 



'/lul. 



2 
-NJRP^, vF^ = --NJR'^K 



(4.47a) 
(4.47b) 



ai 



NJ 



4ao 



3; 

A 



^T>"'^Tj^^ + -T^T^ - ^^2 + ^E - =-A 



1 
16 



2ai Oi 



K, 1 /\/7r,, 
/ ' -)(- ' 



1 21 



+ —YRo 



''TX 



Aaf J 63 6m ai J 12m ai J 



-NJ 






.2^(^Mr^.^gaa)y^y_^ 



(4.47c) 



(4.47d) 



-NJ 



P 



ai 



2ao 

TT -/TTI „|T7\ 



rp{lJ.\cr±rp\iy)_ 



/2 ^_^_ P Ptt _ ^ 

3ao 36ao J 






/2 P, 



"K 



9ao J 



s/ (/.|ra±ga|^>Ny^ 



a-' tctI. 



(4.47e) 



4.3 Massless spin-0 mode 



In the spin-0 mode there is a special "massless" subcase given by 



m = 0, i.e., ai = cq. 
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We take the other parameters in the massless case to be the same as those of the 
massive spin-0~ mode. Therefore we will only mention the parts that differs from 
those in the massive spin-0~ case. The supermomenta and Lorentz rotation parts 
can be further simplified 



n: 







1 %. 



ao 



-. ^ -r ^tlO ,0 J~i^l], 



12 J 



M-L 






n 



l^u 



\T, 






0, 



'TT; 



^'l^r- 



0. 



(4.48) 
(4.49) 
(4.50) 



The non-zero PB's for the PIC's that remain, eqs ( |4.34| , [1.37| ), are {4>Tn7±j 
{(p-p, 4''uB'} ■ Whether the PB's matrix is singular is determined by %/J. Let us first 
suppose the % does not vanish (except perhaps on a set of measure zero). Then the 



class of ' 



>-p±, (Pjlu, 



and 0jjj7_|_ won't change. The SC's are 



X,J[uJ_- 
A.ajn7 - 



J 

T « 



(4.51) 
(4.52) 
(4.53) 



The Lagrange multipliers are 



W^ 



-NJT^", u^"" = 0, 
2 



u' 



s^a„u 



2 2 — 



u 



fii/i. 



P IT 



NJ 



3ai J 



V^ 



<x(A'|ga|i.)y^^^ + r/'^^'^le^l'^^VaT^ - R^'"'^ 



(4.54a) 
(4.54b) 
(4.54c) 
(4.54d) 



We find that x± commutes with 0_l, so its consistency condition leads to, instead 
of determining «-*-, a tertiary constraint (a quite unusual occurrence): 



""n 



«-^7-2;:'«-'=°- 



(4.55) 
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However, (±_ still commutes with (f)± thus «-*" remains undetermined. Then the con- 
sistency condition of C±, i-e., C± ~ 0, and the substitution of the known multipliers 



in ( |4.54a[ -d) lead, in turn to the constraint 

e^ = "RoL ~ 0. (4.56) 



(Because of (|4.49| ) the constraint ^± ~ can be derived directly from the Bianchi 
identity V[aTijc\^ = R[abcf-) Then there is no need to go further, since the constraint 
C_L, along with ^^, contradicts our assumption that % does not vanish. Consequently, 
for the masless spin-0~ case, % must definitely vanish. 

Knowing now that % definitely vanishes, so then does the rhs of ( |4.38| , ^l39| , [4.40| ) , 



which tells us that the massless spin-0 doesn't have any propagating torsion modes 
and is, consequently, just equivalent to GR. 

5 Discussion 

In the PGT, there are forty dynamic variables coming from sixteen tetrad com- 
ponents and twenty-four connection components. The number of total variables 
counts eighty because the same number of canonical momenta accompany the vari- 
ables. Nonetheless, constraints eliminate many unphysical variables. The super- 
Hamiltonian H±, super- momenta TCa, and the Lorentz generators H^u are (generally) 
ten first-class constraints. There are also ten "sure" first-class primary constraints. 
This total of twenty first-class constraints (because of their gauge nature) offset forty 
variables — in general. 

In the dynamic spin-0"'' and massive spin-0~ modes that we have considered here, 
we found that the twenty-three PIC's and eleven SC's are second-class. They thus 
eliminate thirty-four unphysical variables. Therefore, there are only 80 — 40 — 34 = 6 
true physical variables. It simply means that the number of degrees of freedom is 
reduced to three: the scalar/pseudoscalar mode and the two helicity states of the 
usual massless graviton. 
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However, for the massless spin-O" mode, we have shown that the pseudoscalar 
canonical momentum must vanish weakly in order to make the theory self-consistent. 
In that case we have the usual 4 + 4 "sure" first class constraints associated with 
the translational gauge freedom and the 6 "sure" first class Lorentz gauge freedom 
constraints ( |3.2| ). All of the "if" constraints ( |4 .271 - 1^^331 ) then turn out to be first 



class. Their consistency conditions, as well as the consistency condition for ( p.2|) , 
degenerate to the single first class condition % ?« 0. Thus we have a total of 
4 + 4 + 6 + 23+1 = 38 first class constraints leading to 80 — 2 x 38 = 4 physical 
initial values for just the usual 2 degrees of freedom for the graviton. In this special 
case the propagating torsion modes are unphysical, being pure gauge. 

One can determine the positivity of the non-zero coupling parameters from the 
Hamiltonian density ( [4.2a| , [4.2b| , 14.26aD . Since the kinetic energy density must be 



positive definite, we have 



ai>0, — <0, — >0. (5.11 



The results are consistent with restrictions that earlier researchers have proposed. 



e.g., reterencegaaffi. 

It is no wonder that the Lagrange multipliers are identical to the missing veloc- 
ities. Let us examine these in the spin-0"'" mode. The time derivatives of canonical 
variables are given in ( [4.23| ). After reducing the results we truly find 

(5.2a) 
(5.2b) 
(5.2c) 

m±o = %. (5.2d) 

The situation of the massive spin-0~ mode is much the same. On the other hand, the 
expressions for the multipliers, the m's, show that each missing velocity is determined 
by non-velocity terms, hence they do not have a physically independent dynamical 



Tx'^^ = M^^, 


jn^lW] = ^W^ 


Rj_'^ = u^, 


R^_^>'^^^ = u>'^^, 


'R±''^'' = 'u"^'', 


R^^^"'^^=u^"'^, 
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status. This fact reflects the details of the Lagrangian field equations. Due to the 
specific parameter combinations ( ^.1| ) and ( f4.25| ), the coefficients attached to the 
time derivatives of those missing velocities vanish and the related field equations 
become constraints. 

If we look at the multipliers for the two spin-zero modes, u^-^, u'^'^, u^, and u'^'^^ 
essentially have the same expressions, (n'"^ vanishes in the spin-0"*" and massless 
spin-0~ modes because Tpp^ ^ there.) It is clear that, since (f)-p±, 4>-pjj, (p-p, and 

A 

(f)-pi;j_ contain the variables related to spin-one modes, they should not be involved in 
the field equations of spin- zero modes. As to the multipliers u^'^^ and hi"^'^ , which 
are related to the spin-2"'" and spin-2~ modes, in both of the modes they are more 
involved and complicated because the usual massless spin-two graviton is mixed up 
here. 



We note that the constraints ( [4.20|) in the spin-0"'" case deduced from the Lorentz 
rotation parts are intriguing. Since w is generically a function of spatial coordinates, 
the constraint will generally prevent txj from vanishing unless it vanishes globally. 
However w could vanish permanently and become a new "constraint" if its initial 
value is zero. (Assuming that T^ cannot become unbounded.) From ( ^.23dD the 
relation between w and the (affine not Riemannian) scalar curvature is given by 

vj = — R — 3m, (5.3) 

therefore w = indicates that the (affine) scalar curvature is equal to Gum/bQ, a 
constant. The constant forms a barrier for the scalar field to cross. Whether the 
scalar mode evolves on such an affine (anti) de-Sitter spacetime background or is 
frozen depends upon the initial conditions. This forms a new kind of constraint bi- 
furcation. The usual type, which originates from the PB's matrix of the constraints 
being not of constant rank, has been linked to acausal propagating modesi^. (Be- 
cause of the mass terms, as can be seen from ( [4.46 )), the massive spin-0~ mode 
avoids getting involved in this new type of constraint bifurcation.) 
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It is straightforward to work out, from the covariant analysis of the scalar PGT 
modesis, the Lagrangian field equations for the spin-0+, case with cu = 0. They just 
turn out to be those of GR with the addition of a corrected cosmological constant 
depending on A and the value of the constant (affine) scalar curvature, Qnm/h^. All 
coefficients of the torsion components in the field equations vanish. Consequently 
the torsion is really pure gauge in this case. A further very special highly degenerate 
subcase, with w = Q = m, needs no additional discussion; it has just the behavior 
expected in the cjj — > limit. 

This covariant analysis is in accord with the corresponding Hamiltonian analysis 
which shows that the 4 + 4 "sure" translational gauge generators, the 6 "sure" first 
class rotation gauge generators and the 23 primary "if" constraints give rise to just 
one secondary constraint: ro ?^ 0. All the "if" constraints turn out to be first class. 
Thus we again have a total of 38 first class constraints and the usual 2 degrees of 
freedom for the graviton plus purely gauge torsion. 

For solutions with special symmetries, of course there are other possibilities. In 
particular, consider the case of homogeneous cosmologies. If w is dependent on 
time only, the constraint ( [4. 201 ) simplifies, thereby allowing for simpler solutions 
for scalar fields on a Robertson- Walker-like spacetime. But here we are concerning 
ourselves with the dynamic structure of the general theory, not with the pecuharities 
of solutions with special symmetries. 

We remark that the fourth-derivative gravity which people have investigated re- 
cently has similar degeneracy problemEJ. Careful consideration is required to under- 
stand if the relation between its degeneracy and constraint bifurcation phenomenon 
exits, but such a study is beyond the purpose of this paper. 
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6 Conclusion 

In this paper the Dirac-"if" constraint Hamihonian formahsm for the PGT has been 
apphed to the study of the full nonlinear behavior of its spin-zero modes with certain 
specific coupling parameter choices. It is shown that by using the Hamiltonian 
analysis one can clearly identify the degeneracy and the corresponding constraints, 
and the true dynamic degrees of freedom of the PGT as well as the positivity of its 
couphng parameters. 

For the spin-0~ mode, there is a massive, non-ghost pseudoscalar field, in addi- 
tion to the usual spin-two graviton, propagating dynamically. But the corresponding 
massless mode has additional gauge freedom in the torsion modes, yielding a dy- 
namics which is essentially equivalent to that of GR (with presumably undetectable 
purely gauge torsion). Therefore, the magnitude of the mass m will determine the 
viability and detectability of this mode. If |r7i| ^ 0, we might expect to find this 
mode at about the Planck-scale range. On the other hand if \m\ -^ 0, the differences 
between the PGT and GR are too tiny to be detected. In both of these situations we 
can only obtain almost the same effects as GR on the ordinary scale. There is also 
a scalar field propagating in the spin-O"*" mode case. This mode encounters a type 
of constraint bifurcation phenomenon, depending on whether w vanishes, which di- 
vides the phase-space into two subspaces, since vu cannot vanish in general unlesss it 
vanishes globaly because of the Lorentz rotation constraint (|4.2(J|). However, as we 



discussed above, the scalar momentum could be only a function of time which turns 
off the effect of the phenomenon. Then we have a solution which can be applied to 
spatially homogeneous cosmological models. 

The constraint bifurcation phenomenon is generally expected to occur in the full 
nonlinear PGT with higher spin propagating modesEj. Somewhat surprisingly, we 
also see the "constraint bifurcation" phenomenon in these spin-zero cases, albeit in 
a very simple essentially "all or nothing" fashion. 
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This work is an important complement to the earher work on PGT scalar modesl^ 
and is a necessary preliminary to our next step: examining the spin-one and spin- two 
modes and comparing the nonlinear results to the linearized ones of PGT. 
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Appendix: Non-zero Poisson brackets 

In order to classify the SC constraints and derive the Lagrange multipliers, the values 
of all PB's should be calculated. Here the results of the calculations of non-zero PB's 
in the spin-zero modes of the PGT are presented. 

Results in the spin-0^ mode 

The followings just show that all constraints in the spin-0^ are second-class provided 

At 



ixv 



It 

rx,^'^0W}^ 2,5,.,.e^,'^^eVV,— , (A.l 

{\m^ 24^ (A.2) 

{XTWL.n"'^(t^'^±} - 5^.'^na{Tie%)Va—, (A.3) 

{X77i7±,0W}~ -^-^np[arh]{-pe\)Vanf, (A.4) 

{Xjn7±,0W±}~ -^^V^{-pVu)T, (A.5) 

{\wjn7, (p'-^} ~ ^-Yl%njiprjuTe%) - tr^r/^eV)] V„n^ (A.6) 

A:r7T 



u 



rp 



{'x^/u'^y'^}^ S,AVwVvTe%)Va — , (A.7) 
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Vxwi7,¥TfKp} ^ -^IXnjK^^V^pT]^) - Xrt^j^riTj^^T^)]. (A.8) 



Results in the spin-0 mode 



{X±, m'^^0/m±} ~ -^ye^"''^M/Tp^± - 2(5,.,meVVa — , (B.l 

{X±, ^'^0M ^ ^y^M.'^^eVV.^ + 4.'y«^T^x, (B.2 

|X±,0±I~ J m, (B.3 

J ai 

{Xjn7±,w'^^0^±} ~ 5..'-^r7^(jie%)V„— , (B.4 

Ujn7±,07^|~ —-^VplaVrjijie v)^an'', (B.5 

U7n7±,0^±l~ —-^iTw{-piTv)T, (B.6 

{\^, 07±} ^ %^[tw?-eV) - t^?M^^eV)]V,n'', (B.7 

{^X^,w'^0y^ <5..'tw?-eV)V, — , (B.8 
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Table 1: Primary 'if -constraints, critical parameter values and masses 



